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ABSTRACT 

As a continuation of our previous studies in 2007 and 2008, new photometric obser- 
vations of the T Tauri star TW Hya obtained by the MOST satellite and the ASAS 
project over 40 days in 2009 with temporal resolution of 0.2 days are presented. A 
wavelet analysis of the combined MOST-ASAS data provides a rich picture of co- 
herent, intermittent, variable-period oscillations, similarly as discovered in the 2008 
data. The periods (1.3 - 10 days) and systematic period shortening on time scales 
of weeks can be interpreted within the model of magneto-rotationally controlled ac- 
cretion processes in the inner accretion disk around the star. Within this model and 
depending on the assumed visibility of plasma parcels causing the oscillations, the 
observed shortest-period oscillation period may indicate the stellar rotation period of 
1.3 or 2.6 d, synchronized with the disk at 4.5 i?© or 7.1 i?©, respectively. 
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1 INTRODUCTION 

The data gathered by the MOST satellite in 2007 and 2008 
for TW Hya - a T Tauri-type star of spectral type K7Ve 
and V ~11 mag surrounded by an accretio n disk - were a 
subjec t of a Fourier and wavelet analysis bv lRucinski et all 
(2008). The most interesting discovery resulted from the 
wavelet analysis, which revealed oscillatory features system- 
atically shortening their periods on time scales of tenths of 
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days. With typical periods of 2 to 9 days, a preliminary in- 
terpretation was through instability phenomena within the 
innermost accretion disk. In order to check whether such 
variations are a constant feature of the star, TW Hya was 
observed by the MOST satellite for the third (and most 
likely last) time in 2009 (Section Additionally, simul- 
taneous ASAS observations were obtained as close to the 
MOST observation in time as possible. In this paper we 
present an analysis of the new photometric data (Section [3j 
and propose a more physical interpretation of the features 
discovered in the wavelet spectrum (Section 2]). 



2 OBSERVATIONS AND DATA REDUCTION 

The optical system of the MOST satellite consists of a 
Rumak-Maksutov f/6 15 cm reflecting telescope. The cus- 
tom broad-band filter covers the spectral range of 380 - 
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Figure 1. The 2009 light curve of TW Hya in flux units, scaled (separately for each source) to unity at the mean brightness level. 
The original MOST data are presented as crosses and their respective median-orbital values are shown as black squares. The ASAS 
observations obtained at hjd = 4894.5 — 4903.9 in the Johnson V filter with 10 minute temporal resolution, expressed in flux units 
assuming mean magnitude V = 11.065, are shown as small dots; their mean values, calculated for intervals of 0.2114 day, corresponding 
to the MOST data collection rate of 3 satellite orbits, are shown as circles. Note the good match between the MOST and ASAS data 
sets where they overlap during the first nine days of the MOST observations. 
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Figure 2. The frequency spectrum of TW Hya obtained from 
the combined, time-averaged 2009 MOST (185 points) and ASAS 

V filter (23 points) data. No statistically significant peaks are 
present above the frequency of 1 c/d. The dotted line shows the 
arbitrarily scaled dependence, while the crosses represent 
errors of amplitudes o,(/). 

700 nm with effective wavelength falling close to Johnson's 

V band. T he pre-launch charact eristics of the mission are de- 
scribed by Walker et al. and th e initial post-launch 
performance bv lMatthews et al.l (|2004l ). 

TW Hya was observed between March 3 and April 12, 
2009, corresponding to the hjd range of 4894.53 - 4934.83 B 
The data were collected during the low stray-light segment, 
lasting typically 30 mins, of (almost) every third satellite 
orbit over 40.3 days. The direc t-imaging mode of the satellite 
was used (| Walker et al.ll2003l ) with individual exposures of 
30 sec. 

The data were reduced in the same way as described by 
1 Throughout this paper: hjd = HJD - 2, 450, 000 



ISiwak et all (|2010l ). Small linear correlations between the 
stellar flux and the sky background level within individual 
orbits (most probably caused by a small photometric non- 
linearity of the electronic system) were also removed. In the 
last step the data were de-trended by means of a low-order 
polynomial fitted to the light curve of a constant star, si- 
multaneously observed with the main target. As a result, 
we obtained a very good quality light curve with median 
error of a single data point of about 0.009 mag. The system- 
atic errors related with the uncertainty of the nonlinearity 
corrections do not exceed 0.005 mag and are inconsequential 
for our investigation. 

In parallel with the MOST observations, TW Hya was 
also monitored b y the All Sky Automated Survey (ASAS) 
|Poimanskilll997l ) in Johnson V filter with temporal sam- 
pling of about 10 min. Due to deterioration of the CCD chip 
used for I/-band observations over the last years, the individ- 
ual data points have rather low accuracy with a median error 
of 0.045 mag. The mean observed magnitude, V = 11.065, 
was used for conversion of the ASAS data to norm alized 
flux units. Similarly as in 2008 (|Rucinski et al.ll2008h . the 
2009 data show good correspondence between MOST and 
the ASAS V-band. For comparison with the MOST obser- 
vations, the ASAS l/-band data points were averaged using 
the sampling rate of the MOST data of 0.21142 d. 

In Figure [T] we present the original and the averaged 
data points of TW Hya obtained by the MOST satellite and 
ASAS in Johnson V filter in normalized flux units. The me- 
dian error of the time-averaged MOST data in flux units 
is equal to 0.0066. The time averaged ASAS F-band data 
have an order of magnitude worse accuracy, of about 0.041 
flux units. For both datasets, the errors include the gen- 
uine, time-averaged, short-timescale physical variability of 
TW Hya observed within individual MOST orbits; thus, the 
rapid variability of the star was lost in the averaging pro- 
cess, simil arly as in the previou s analysis of the 2007 and 
2008 data (|Rucinski et alj|2008h . 
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Figure 3. The wavelet transform of the 2009 TW Hya MOST data supplemented by the ASAS data, as explained in the text. The 
amplitude power of the transform is expressed by the grey scale intensity plus contours with arbitrary units. Edge effects are present 
outside the white broken lines; they do not affect our conclusions. At the bottom, the combined MOST and ASAS-V light curve, re- 
sampled into a uniformly distributed time-grid with 0.2114 d spacing, is shown; at the right a comparison is given of the restored (via 
horizontal averaging) frequency spectrum (the continuous line) with the power spectrum corresponding to the amplitude spectrum in 
Fig. [2] (the dotted line). 



3 ANALYSIS OF MOST 2009 DATA 
3.1 Fourier transform 

We performed a Fou rier analysis of the 2009 data as in 
iRucinski et al.l l|2008h (see Section 4.2 of that paper for de- 
tails) . The frequency spectrum was obtained using the light 
curve formed from the averaged MOST and ASAS V-band 
data, represented by squares and circles, respectively, in Fig- 
ured] The 0.2114 d sampling rate (3x the MOST orbit dura- 
tion) precludes analysis of frequencies higher than 2.35 c/d. 
In Figure [2] we show only the frequency range below 1 c/d 
because there are no stati stically signific a nt pe aks at higher 
frequencies. We note that lGiinther et ahl (|2010h did not find 
any sub-minute periodicities in their high-speed time-series 
photometry of TW Hya obtained with the SALT telescope. 

Because the temporal resolution of the 2009 data was 
3x poorer than that of the 2008 data, the amplitude errors 
are now larger than before, with median value of the error 
0.012 of the mean flux. As observed before, the peak ampli- 
tudes a(f) seems to scale as flickering noise (a(f) ~ 1/VDi 
but due to the large r errors, this tendenc y is not as obvious 
as in the 2008 data jRucinski et alj|2008h . 



3.2 Wavelet analysis 

The original time-distribution of averaged data points dif- 
fers slightly from the uniformity required by the wavelet 
technique. Therefore, we interpolated the data into a uni- 
form time grid using splines. Except for the orbital-averaged 
MOST data points, we also included the first five aver- 



age ASAS points obtained in V-band before and during 
the gap in the MOST data which occurred within hjd ~ 
4895 — 4896.5. The remaining ASAS points were not used 
because the ASAS data were significantly less accurate and 
were used here only to assure uniformity of the main MOST 
sequence for the wavelet analysis. 

Figure [3] shows the result of the wavelet transform 
of 191 data points uniformly distributed in time at the 
spacing of 0.2114 d. As in the case of the 2008 data, the 
Morlet-6 wavelet provided the best match between the time- 
integrated power spectrum and original frequency spectrum, 
as is visible in the right panel of Fig. [3] We stress that only 
for the Morlet-6 wavelet do we obtai n such a nearly-perfect 
match of both plots; as described in IRucinski et al.l (|2008l ) 
Morlet transforms of other orders result in systematic dif- 
ferences in the period scale. 



4 DISCUSSION 

The most important result of this study is that the gen- 
eral behaviour of the Morlet wavelet transform in the time 
- period plane i s simi lar to that obtained and described by 
IRucinski et al.l l|2008h for the 2008 observing season. The 
systematically migrating oscillatory features that we saw in 
2008 were not an unusual phenomenon limited to this par- 
ticular point in time; apparently, TW Hya shows this phe- 
nomenon all or most of the time. The new 2009 data not only 
confirm the direct visibility of the features drifting to shorter 
periods, but this time we could identify as many as four of 
them. Their approximate moments of appearance and dis- 
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Table 1. Numerical values extracted from the wavelet spectrum 
and related physical parameters computed assuming that the fea- 
tures observed in the spectrum come from a Keplerian disk around 
the central star. 

The features marked with an asterisk * arc for the data in 
iRucinski etafl (I2008T) . 
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-0.065(1) 


5.2 


8.2 


3 - start 


4902 


7.2 




14.1 


22.3 


3 - end 


4931 
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6.0 




12.5 


19.8 


5* - end 


4553 


3.0 


-0.143(2) 


7.8 


12.5 


6* - start 


4556 


7.5 




14.5 


22.9 




4900 4910 4920 4930 

HJD - 2,450,000 

Figure 4. The linear fits (thick black lines) to the local maxima 
of features 1-3 visible in Fig. [3] Similar fits were made for 
the features observed in 2008 of which only feature no. 5* could 
provide sufficiently well defined results. 



appearance (in hjd) , the periods P observed at the beginning 
and at the end of each oscillation packet and the respective 
Keplerian radii m calculated assuming M* =0.72 Mq and 
R, — 1 i?Q are presented in Table [1] they are ordered in 
their start periods. The data in Table [1] for features no. 1-3 
were obtained by tracing local maxim a along the patterns in 
Figure [3] similarly we used figure 8 of IRucinski et all (|2008l ) 
to trace feature no. 5* identified in the 2008 data. We did 
not investigate features nos. 4 and 6*: the former had a long 
period of almost 10 d, so that the 2009 run was too short to 
fully follow its evolution while the latter feature (from the 
2008 data) occurred too close to the end of that observing 
run. 

Within the accuracy of locating the local peaks, only the 
lowest-order fits by straight lines (see Figure 3} gave a suf- 
ficiently significant description of their temporal evolution 
described by their time derivatives P, as given in Table Q] 
The fits were weighted by their local peak intensities. We 
note that the rate of the period change appears to be faster 
for longer periods. 

In the now most popular interpretatio n, originally 
develo pe d for accreting neutro n stars by I Ghosh et al.l 
1 19771) ; iGhosh fc Lambl (|l979al lbh and fu rther applied 



to classical T Tau- ty pe sta rs (CTTS) 
Cameron & Campbell! Il993t iHartmann et al 



(iKonigll 



Shu et al 



19941 ; 



Muzerolle et all 199 



2001] ) 



1991 



1994 



magnetic 

loops originating on the stellar surface connect to the 
accretion disk outside the co-rotation radius to avoid 
angular momentum transfer from the disk. The model of 
magnetically controlled accreti on explains observa tions 
of many CTTS, e.g. AA Tau (|Bouvier et al.l 120071 ) and 
BP Tau (jDonati et al.ll2008r i. as well as some young brown 
dwarfs, e.g. 2MASSW J1207 334-393254 in the TW Hya 

association (IScholz et al.ll2005h. 

As proposed before l|Rucinski et al.ll2008h . we think that 
the observed light oscillations are produced in one or sev- 
eral regions in the inner accretion disk due to instabili- 
ties in the disk plasma. They may be caused by magnetic 
field threading of the disk near its inner edge via Kelvin- 



Helmholtz instab ility or by combined t urbulent diffusion 
and reconnection. IGhosh fc Lambl (|l979al lbh found that such 
an interaction would produce a broad transition zone be- 
tween the unperturbed Keplerian disk and the co-rotating 
magnetosphere outside which considerable departures from 
Keplerian motion are expected. Their theoretical expecta- 
tion for the o uter t ransition zone seems to be confirmed by 
ICalvet et all (|2002r ) (their Sec. 4.3), who found that at a 
distance of 1 AU, gas in the TW Hya disk probably drifts 
slowly inward nearly following Keplerian orbits. In this pic- 
ture, the oscillations observed by MOST may be explained 
by orbital motion of optically thick condensations produced 
by magnetic loop - disk interactions in the inner transition 
zone within more spatially confined regions, bright enough 
to produce the observed, large amplitudes of light variations. 
Once initiated at respective Keplerian radii a\ (see Table [TJ, 
they would then rapidly propagate toward the co-rotation 
and truncation radii, a property which could explain the 
high absolute values of P. 

Although there may be doubts if the stellar rotation 
- disk revolution synchronization still operates in TW Hya 
(see below), one could this way estimate the period of ro- 
tation of the star. The shortest period observed in the 
MOST data (the feature no.l at the end of its evolution at 
hjd = 4913) is 1.3 ± 0.1 d; admittedly, this period is poorly 
defined at the 2.2er significance level. This value is almost 
identical to the shortest period amo ng many periods previ - 
ously claimed for TW Hya of 1.28 d (frlerbst fe KoretJll988r i 
(see sect. 2 of IRucinski et al.l i|2008T ) for a summary of this 
topic). Simi lar values (1.4±0.1 d and 1 .55zfc0.10 d) have been 
obtained bv lAlencar fc Batalhal (|2002T ) in an analysis of veil- 
ing measurements, but the authors disregarded these values 
paying more attention to other detect ed periodicities, i.e . 
2.9 and 4.5 d. However, as obtained bv lWilner et all (|2000l ) 
from 7 mm VLA interferometric observations, the outer ac- 
cretion disk is oriented almost face-on. This means - for the 
inner accretion disk being co-planar with its outer parts - 
that each hypothetical condensation causing observed light 
oscillations is seen twice during each orbital evolution. Then, 
the observed periods given in Table [1] would in fact repre- 
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sent half-values of the real orbital periods. Accordingly, we 
re-calculated the Keplerian distances assuming doubled val- 
ues of the observed periods, which are given as 02 in the last 
column of Table [1] In this scenario, the value of 2.6 ± 0.2 d 
rather than 1.3 ± 0.1 d, would correspond to the pure rota- 
tional period of the central star itself. 

Several semi-direct esti mates of the inner d i sk rad ius are 
available from the literature: IWeinberger et all ((2002), using 
spatially resolved HST and Keck infrared images and spec- 
tra, modeled assuming different dust grain sizes, obtaine d 
10.7 ±2.1 Rq for the inner disk radius. lEisner et all l|2006h . 
using spatially resolved near isT-band interferometric Keck 
observations and a model utilizing sub micron-sized dust 
grains , obtained 12.9±2.1 Rq. Contrary to lWeinberger et al.l 
(2002), their second model with large dust grains gave a 
smaller radius of 4. 3 ±2.1 Rq, but much poorer fit to all the 
data. We note that only if we assume the doubled orbital pe- 
riod do we obtain an inner disk radius of 7.1 ±0.4 Rq, which 
would be commensurate with these estimates (Table [TJ. 

In spite of its conceptual appeal, the coupling between 
stellar and disk rotation offered by the magnetically con- 
trolled accretion model may no longer be valid for TW Hya. 
Although TW Hya possesses properties of a young (1-2 Myr 
old) CTTS, it is part of an evolved (~ 12 ±8 Myr old) TWA 
association, as obtain ed from lithium depletion estimates of 
iMentuch et al.l (|2008l ). The accretion disk itself also shows 
indications of bei ng in advanced stages of dust evolution: as 
previo usly noted. [Weinberger et alj (|2002l ) and lCalvet et al.1 
( 2002) needed larg e-size dust grains to fit the observed spec- 
tral energy distribution; the latter authors discovered also a 
developing gap in the disk at the distance of 4 AU from the 
central sta r. From an analysis of sev eral young stellar as- 
sociations, IJavawardhana et al.l (|2006l ) found that accretion 
disks in T Tau-type stars appear to cease before the age of 
about 10 Myr: among 32 members of the TWA association, 
only TW Hya and TWA 3A show signs of ongoing accretion 
and seem to be slow rotators. We note that the small value of 
v sini = 10.6 km/s determined by the authors for TW Hya 
should not be taken to indicate a long rotational period of 
the star, but is rat her a result of the low inclination angle 
of its rotation axis (|Wilner et all 120001 ). 

Because of the p ole-on geometry of TW Hya, 
lAlencar fc Batalhal (|2002l ) conclude that it is rather unlikely 
that the observed light variations are caused by rotation of 
hot spots on the star photosphere, especially in the case 
when magnetic poles seems to be aligned with the stellar 
rotation axis (see Section 4 of their paper for detailed expla- 
nation) . These considerations may also apply to the theoret- 
ical prediction of regularities in the geometrical spot distri- 
bution as well as quasi-periodic ligh t oscillations caused by 
hot s p ots on the stellar photosphere l|Kulkarni fc Romanoval 
120091 : iRomanova fc Kulkarnl 120091 ). Although the oscilla- 
tions may have similar characteristics to light variations ob- 
served by MOST, a long rotational period of TW Hya (of 
order 10 days) would be required to be explained by such 
a mechanism. We note that neither the 2008 nor the 2009 
MOST observations revealed any common for both runs, 
long-period (8-15 d) wavelet-transform features which could 
be attributed to rotational period of the star itself. Although 
this result may be caused by the limitation of a single MOST 
monitoring run to about 45 day, there are no indications of 



such a long rotationa l period in the Fourie r spectrum of the 
long series of AS AS (|Rucinski et aljfeoosl ) observations. 



5 CONCLUSIONS 

The wavelet transform analysis of the 2009 data obtained 
by the MOST satel lite and the ASAS pr oject confirms the 
essential findings of lRucinski et al.l (|2008l ): The general red- 
noise nature of the Fourier spectrum of TW Hya is rather 
complex in that it corresponds to several oscillation packets 
with their frequencies increasing in time and with larger am- 
plitudes for lower frequencies. The respective orbital radii 
can be calculated assuming that periods of individual os- 
cillation packets are related to Keplerian motion of op- 
tically thick plasma condensations in the inner accretion 
disk. The observed periods are consistent with their ori- 
gins in the inner accretion disk, several stellar radii out- 
side of the co-rotation and truncation radii. Noting the 
pole-on geometry of the accretion disk for TW Hya, con- 
sideration should be given to the possibility that the ob- 
served periods represent the respective halves of the true 
orbital periods. Assuming that we detected the shortest pos- 
sible period in the Morlet wavelet transform of 1.3 d and 
that the inner boundary of the disk co-rotates with the 
star, the co-rotation radius for TW Hya would be about 
4.5 ± 0.2 Rq or 7.1 ± 0.4 Rq depending on the assumed 
period (1.3 ± 0.1 or 2.6 ± 0.2 d, respectively); only for the 
longer pe riod is the inner disk siz e in accordance with th e 
results of IWeinberger et al l (|2002l ) and lEisner et al.l [|2006h . 
obtained from spatially-resolved infrared images and inter- 
ferometric observations. Determination of the rotational pe- 
riod of the star is crucially important both for the inner ac- 
cretion disk consideratio ns as well as for mechanisms i nves- 
tigated theoretically by iKulkarni fc Romanoval (|2009l ) and 
IRomanova fc Kulkar ni (2009) which could explain light vari- 
ations through (possibly time-evolving) geometrical distri- 
bution of hot spots and/or quasi-periodic oscillations in this 
region. 
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